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1-001 ° ° pili green 1o8d 1i3BTy EAYANTS eaesoreeoucu.mars1 2% cy 


OWDLU = D_floating ** unsigned posener pover routine 
V/ ; File: OTSPOWDLU.MAR Edit: SBL1001 


PARA ARARAAAAAAALAAALALL ALLELE ESET RESTS TREE EE EES ETE E CECT PCP ce cece ccc aca? | 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
ONLY He ACCORDANCE Me THE TERMS OF SUCH obtcane’ wee arch 


IN 0 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 


® fe 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
CORPORAT I Ow NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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+ 

FACILITY: Language support procedures, Mathematics Division 

ABSTRACT: 


This module contains OTS$POWDLU, a procedure which raises a 
D_floating value to an unsigned integer power. 


ENVIRONMENT: Runs at any access mode, AST Reentrant 
AUTHOR: Steven B. Lionel, CREATION DATE: 22-JUL-1981 
MODIFIED BY: 

_1-001 - Original. SBL 22-JUL-1981 
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1-001 DECLARATIONS . wis bree 1 Oba 94:38:44 MTHRTL.SRCJOTSPOWDLU.MAR; 1 ° 3) 
p ¢$ «SBTTL DECLARATIONS 
rk ; LIBRARY MACRO CALLS: 
0 SSSDEF 3; System error codes 
1 SCHF DEF ; Define condition handler symbols. 
¢ SSFDEF 3; Define stack frame symbols. 
4 SPSLDEF ; Define program status longword 
$3 2 3 symbols. 
0000 >% 
44 3? ; EXTERNAL DECLARATIONS: 
$008 8 DSABL GBL : Force all external symbols to be declared 
0000 59 <EXTRN MTHSK_UNDEXP 3; Undefined exponentiation 
0000 60 -EXTRN MTHSK_FLOOVEMAT ; Floating overflow in Math Library 
0000 61 ~EXTRN MTHSK_FLOUNDMAT ; Floating underflow in Math Library 
0000 66 -EXTRN MTHSSSIGNAL ; Math error routine 
0000 635 ; 
0000 64 ; MACROS 
0000 65 ; 
0000 66 ; NONE 
0000 67 ; 
0000 68 ; EQUATED SYMBOLS: 
0000 69 ; 
0000 70 ; NONE 
0000 713 
0000 ig: ; OWN STORAGE: 
0000 is 3 
0000 74 ; NONE 
0000 9 § 
B08 6 3 PSECT DECLARATIONS: 
00000000 78 -PSECT _OTSSCODE PIC, USR, CON, REL, LCL. SHR, - 
s'4'4 EXE, RD, NOWRT, LONG 


No 


Mm 14 


; CALLING SEQUENCE: 
result.wd.v = OTSSPOWDLU (base.rd.v, exponent.rlu.v) 
; FORMAL PARAMETERS: 


SOOoooooooooooo 
COOoooooooooooooo 


00000004 base = 4 floating base. Note that this is passed b 


l 
mmediate value in two argument List positions, 
c 
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1 1 OTSSPOWDLU : D_floating ** unt igned inte 6-SEP-1984 94:38:49 UMTMRTL. SREIOTSPOWDLU.MAR: 1 a 3) 
p § : -SBTTL OTSSPOWDLU - D_floating ** unsigned integer power routine 
4 ; FUNCTIONAL DESCRIPTION: 
§ : Raise a D ay tte, value to an unsigned integer power giving 
o ; a D_floating result. 
3 ° 3 The result is given by: 
0 1: base exponent result 
8 4 ; any >0 product (base * 2**i) where i is each 
4 ¢ : non-zero bit position in exponent 
$000 g 4 + Undetined tati 
$ = = ndefined exponen on 
0000 «98: <6 ef 1-0 vee sop ne 
8 0 9; 
0 93 
0000 13 
0000 § 3 
0000 é 
0000 4; 
23 
6 3 
7 
8 
9 


hich is a violation of the calling standard. 
his is allowed for language support procedures. 


0000000C exponent = 12 ; Unsigned longword integer exponent 


IMPLICIT INPUTS: 
NONE 
IMPLICIT OUTPUTS: 
NONE 
ROUTINE VALUE: 
The base raised to the exponent's power. 
SIDE EFFECTS: 
May signal: 
MTH$_FLOOVEMAT = floating overflow in Math Library 
MTHS_FLOUNDMAT = floating underflow in Math Library, only 
the caller has FU enabled. 


MTHS_UNDEXP - undefined exponentiation, if both base and 
exponent are zero. 
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OTSSPOWDLU - Of 
1-001 OTSSPOWDLU = D_floating ** unsigned inte 6-SEP-19 


0 V04-00 
S 


r a 4 oT: 
MTHRTL. SRCJOT = 2=( 


Pp 
POWDLU.MAR; 1 (3) 


001¢ ; 133 ENTRY OTSSPOWDLU, “M<R2,R3,R4> ; Entry point 
6D 64°AF 9 141 OVAB mis HANDLER, (FP) ; Translate exceptions to MTH$ errors 
9940 8F 8 6 146 BISPSW #PSLSA_FU ; Enable floating underflow detection 
0 08 pA 14 OvD . ; Ad = initial result 
2¢ 04 AC 7 D 144 MOVD base(AP), R2 ; R2/R3 = base 
4 OC aC OD 11. 145 MOVL exponent (AP), R4 3; R4 = exponent 
93 1 8 13 146 BNEQ BEGIN 3; Skip if not zero 
3 7 01 147 TSTD ; Is base also zero? 
2 1 Bolg 133 BEQL UNDEF INED : jf 80. error. Otherwise, answer is 
04 bpi8 139 RET ; With result of 1.0 
1 Way 
001C 13 3+ 
Q01C 153 ; Do the first iteration here so that we can clear the high-order exponent 
itt: 13 ; a Afterwards, that bit will be zero so we don’t have to worry about 
001C 156 ;:- 
001C 157 
001C 158 BEGIN 
0354 00 €5 OO1C 159 BBCC #0, R4, 10$ : Skip initial multiply 
50 52 70 0020 160 MOVD R2 ; New partial result 
—> Ss F&F 9C 0023 161 108: ROTL #-1, RG, RG ; Get new low exponent bit and clear 
0028 196 : high bit 
10 13 0028 16 BEQL DONE : Is that all? If so, we're done. 
002A 164 
002A 165 ;+ ; 
002A 166 ; Each time we get here we know that there's at least one more exponent bit 
002A 167 ; left, so square our current power of the base. 
002A 168 ;- 
002A 169 
002A 170 LOOP: 

52 52 64 QO2A 171 MULD2 R2, R2 ; Square base. This may overflow or 
002D Ze ; underflow. In either case, the 
b055 17 : final result would also overflow 
0020 174 : or underflow. Our exception 
002D 175 : handler will catch these cases. 

03 54 =€9 BASS 176 BLBC R4, NEXT_BIT : See if this exponent bit is on. 
0030 177 : If so, multiply base by product. 

50 52 64 0030 178 MULD2 R2, RO 3; This too ey overflow or underflow, 
0033 +179 :; which will be caught by the handler. 
Oos3 180 NEXT_BIT: 

54 54 FF 8F 78 3 181 ASHL #-1, R4, R4 3; Get next bit of exponent. 
FO 12 $38 186 BNEQ LOOP 3; Keep going if not all done. 
Beh 184 DONE: 
04 QO3A 185 RET ; Return with answer in RO/R1. 
bots 186 
038 187 ;+ . 
ote ' 3 3 Undefined exponentation error - 0**0 
$3 190 
038 #191 UNDEFINED: ; ; i 
00'8F 9A B 136 MOVZBL #MTHSK rs -(SP) ; Indicate undefined exponentiation. 
50 01 OF 79 F 19 ASHQ #15, #7, R ; RO/R1 = reserved floating operand 
15 43 19% BRB SIGNAL_ERROR 
4 195 
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1-001 OTSSPOWDLU D_ floating ee unt gned inte o- SEE} 984 1338343 MTMRTL.SRCJOTSPOWDLU.MAR;1 © 3d 
4 196 ;+ 
r 199 ; EXC_HANDLER unwinds here if we got a floating overflow. 
23 198 | 
4 00 OVERFLOW: 
7E OO°8F 9A Se 01 MOVZBL #MTHSK_FLOOVEMAT, =(SP) ; Floating overflow in Math Library 
~~ Oo Fe FF 49 8 ASHQ . R ; RO/R1 = reserved floating operand 
0B 11 0040 20 BRB SIGNAL_ERROR 
O4F 04 
O4F 09 ;+ 
O4F $ ; EXC HANDLER unwinds here if we got a floating underflow. 
Oper 8 ; If our caller enabled FU, then Signal FLOUNDMAT, with a result of zero. 
O4F 8 ; Otherwise quietly return zero. 
Beer 09 ;- 
O4F 10 
OO04F 11 UNDERFLOW: 
50 7C QO4F ¢ LRD RO ; Initial result 
E404 AD 06 €1 0051 1 BBC #PSLSV_FU, SFSW_SAVE_PSW(FP), DONE ; Just return if 
0056 14 ; caller disabled underflow. 
7E O00O°8F 9A 0056 15 MOVZBL #MTHSK_FLOUNDMAT, -(SP) ; Floating underflow in Math Library 
005A 18 BRA SIGNAL ERROR 
005A 1 
OOSA 18 3+ 
OOSA 19 ; Signal a MTH$ error whose code has already been pushed on the stack. 
005A 220 ; RO/R1 has either a reserved operand or a zero, depending on the 
OOSA $ 1 ; exception type. 
OOSA ¢ i- 
OOSA 
OOSA 4 SIGNAL_ERROR: 
6D D4 OOSA 5 CLRL (FP) ; Cancel our condition handler | 
00000000'GF 01 FB 005 § 6 CALLS #1, G*MTHSSSIGNAL 3; convert to 32-bit condition code 
006 27 ; and signal a MTH$ error. Omit | 
0063 228 ; second argument to show that this 
006 4 ; is a CALL entry point. : 
04 006 § 0 RET ; Return to caller with result in 
0064 31 ; RO/R1 
0064 232 | 
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“M<> 


#SS$_RESIGNAL, RO 
CHFS$L_SIG_ARGS(RO), R1 
B*OVERFLO@, -4(R0)ER1) 
CONTINUE 
CHF$L_SIG_ARGS(RO), R1 
B*UNDERFLOW, -4(ROSCR1J 
CONTINUE 


#SS$_CONTINUE, RO 


floating ** unsi 
POWDLU = D_ 
64 4; 
Se Sa: EXC 
; while per 
64 3 5 
Bat eles: atten 
é erwise resignal. 
064 40 ;:- ° 
064 4) 
064 t EXC_HANDLER: 
064 4 WORD 
066 44 
s44 45 MOVL 
A 6 TSTL 
06D 4 BNEQ 
Boer 48 MOVL 
0 ; 49 MOVL 
07 50 CMPW 
07C 51 BEQL 
443 26 CMPW 
008 5 BEQL 
0085 54 CMPW 
008A 55 BEQL 
008C 2$ CMPW 
0091 5 BEQL 
0093 58 RESIGNAL: 
0093 59 MOVZWL 
0098 60 RET 
0099 61 
0099 6¢ GOTO_OVERFLOW: 
0099 6 MOVL 
009C 64 MOVAL 
QOA 65 
O0A 66 BRB 
QOAS 67 
QOA4 68 GOTO_UNDERFLOW: 
OOA4 69 MOVL 
Q0A7 70 MOVAL 
OOAD 71 
OOAD i : BRB 
OOAD 7 
QOAD 74 CONTINUE: 
OOAD 75 MOVZWL 
0080 76 T 
081 77 
081 78 
0081 79 ~ END 


Th 28:41 


+ 

EXC eo = - This condition handler gets control if we got an exception 
orming the exponentiation. If it is a floatin 
then continue execution at OVERFLOW, to cause MTH$_FLOOVEMAT to be 


Similarly, continue at UNDERFLOW if an underflow was seen. 


; Get number of si 


AX/VMS Macro V04-00 
MTHRTL.SRCJOTSPOWDLU.MAR; 1 


overflow, 


Entry point 


Get mechanism arguments List 
At depth zero? 

If not, resignal 

Get signal arguments List 
Get signal name 

Overflow fault? 

If so, continue at OVERFLOW 
Overflow trap? 

If so, continue at OVERFLOW 
Underflow fault? 

If so, continue at UNDERFLOW 
Underflow trap? 

If so, continue at UNDERFLOW 


Resignal exception 
Return to VMS exception dispatcher 


Get number of signal arguments 
Move address of OVERFLOW routine 
to PC in signal argument List 
Continue execution 


anal arguments 
Move address of UNDERFLOW routine 
to PC in signal argument List 
Continue execution 


Continue at PC in signet List 
Return to VMS exception dispatcher 


End of module OTSSPOWDLU 
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Symbol table 6-SEP=19 1:28:11 CMTHRTL. SRE ¢i0 OTSPOWDLU.MAR; 1 (4) 


BASE = 00000004 
BEG it R 02 
CHF SL_MCHARGLST = 
CHFSL-MCH_DEPTH = 
CHF SL_SIGARGLST = 4 
CHFSL_SIG_ARGS = 0 
CHF SL~S1G-NAME = 4 
CONT IRUE 83 AD R 0 
DONE 0 3A R 8 | 
EXC_HANDLER $ 4R 
EXPONENT = 0 4 
GOTO_OVERFLOW ° 3 R 8 
GOTO™ .UNDERF LOW AS R | 
L 000002A R 8 
MTHSS$SIGNAL reereree x 
MTHSK_FLOOVEMAT reeeeeee x 00 
MTHSK_FLOUNDMAT eeeneeee XxX 00 
MTHSK_ UNDE XP eereereer KX 00 
IT 0000033 3 
OTSSPOWOLU 0000000 RG 0 
ERFLOW 0000004 02 
PSLSM_FU = 4448 eh 
PSL$V_FU = 0000 $06 
RE AL 0000009 02 
SF$W_SAVE_PSW = 00000004 
SIGNAL_ERROR 0000005A 02 
SS$_CORT INUE = 00000001 
SS$_FLTOVE = 0000048C 
SSS_FLTOVF _F = 00000484 
SS$_FLYUND = 0000049C¢ 
SS$_FLTUND_F = 000004C4 | 
SS$_RESIGNAL = 00000918 
UNDEF INED 00000038 R 02 
UNDERF LOW OOO0004F R 02 
bun wer onrwe ner woranne + 
: Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
. . 00000000 ( 99 ; 9. ) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
SABSS$ 00000000 ( 0:) 4°} NOPIC USR CON ABS LCL NOSHR EXE RD WRT NOVEC BYTE 
_OTSSCODE 00000081 ( 177.) 3 2. PIC USR CON’ REL LCL SHR EXE RD NOWRT NOVEC LONG | 
Seeman re sens an eanmreaese seo on | 
H | Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 29 90 :00:00.09 9:00:01 0) 
Command processing ig : 8: Q. :00:03.14 
Pass 1 21 :00: 9 : he 82 
Symbol table sort 0 :00:00. :00: 1-6 
Pass 61 :00:01.1 :00: 3° 4 
Symbol table output 6 :00:00.05 :00:00.13 


sche pene ee } 
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VAX=-11 Macro Run Statistics 6-SEP-1984 11:28:11 CMTHRTL.SRCJOTSPOWDLU.MAR; 1 (4) | 


uate = BRE 


ine york ing seh Limit was 1200 pages. 
26577 bytes_(52 pages) of virtual memory were used to buffer the intermediate code. 
There were 30 pages of symbol table space allocated to hold 508 non-local - 1 local symbols. 
source Lines were read in Pass 1, pred (ne 13 object records in Pass 2. 
n 


1 pages of virtual memory were used to define 10 macros. 


$we reece ree te ecs ee nese ese meee > 


! Macro Library statistics ! 


| 
teocecococos eooeeoeceocccs - | 


Macro Library name Macros defined 


562 GETS were required to define 7 macros. 
There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESSION/D1I SABLE=(GLOBAL , TRACEBAcK) /LIS=LIS$:OTSPOWDLU/OBJ=OBU$:OTSPOWDLU MSRC$:OTSPOWDLU/UPDATE=(ENH$:OTSPOWDLU) 


244 AH BT13A-SE : ENT CORPORATION 
VAX/VMS V4.0 ND PROPRIETARY 


